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Abstract 

The focus is on a model reduction framework for parameterized elliptic eigen¬ 
value problems by a reduced basis method. In contrast to the standard single output 
case, one is interested in approximating several outputs simultaneously, namely a 
certain number of the smallest eigenvalues. For a fast and reliable evaluation of 
these input-output relations, we analyze a posteriori error estimators for eigenval¬ 
ues. Moreover, we present different greedy strategies and study systematically 
their performance. Special attention needs to be paid to multiple eigenvalues 
whose appearance is parameter-dependent. Our methods are of particular inter¬ 
est for applications in vibro-acoustics. 


1 Introduction 

For the fast and reliable evaluation of input-output relations for parameterized partial 
differential equations (/iPDEs), reduced basis methods have been developed over the 
last decade; see, e. g., [34, 36] or [33, Chapter 19] for comprehensive reviews, with the 
first reduced basis problem being investigated in the 1980’s [28]. The methodology has 
been applied successfully to many different problem classes both in the real-time and 
the many-query context. These problem classes include among others finite element 
discretizations of elliptic equations [36], parabolic equations [15, 35, 39] and hyper¬ 
bolic equations [7, 16]. Furthermore the reduced basis method has been extended to 
Stokes problems [20, 23, 37, 38] as well as to variational inequalities [13, 17] with a 
time-space formulation of the problem and corresponding analysis. It was also applied 
to Stochastic Processes in [11,41] and to a finite volume scheme of a parameterized and 
highly nonlinear convection-diffusion problem with discontinuous solutions in [10]. 
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A posteriori error estimators w. r. t. parameter variations generally facilitate the con¬ 
struction of reduced basis spaces by greedy algorithms as well as the certification of 
the outputs of the reduced models. Different greedy methods for reduced basis and 
error estimators have been introduced in [26, 27, 40], also a greedy method for eigen¬ 
values is introduced in [6], and the convergence of greedy methods has been analyzed 
in [4, 5, 8]. Another way to construct a reduced space is the proper orthogonal decom¬ 
position (POD) method, as discussed, e. g., in [21, 33], 

The problem class of parameterized elliptic eigenvalue problems (/iEVPs) is highly 
important but up to now only marginally investigated in the context of reduced basis 
methods. The first approach [24] from the year 2000, which is based on [25] among 
others, is restricted to the special case of an estimator for the first eigenvalue. In the 
following publications [30-32], the method from [24] is developed further to include 
several eigenvalues. 

However, both the analysis and the algorithms do not cover the case of multiple 
eigenvalues. Quite often, the “vectorial approach”, i. e., the treatment of the eigen¬ 
vectors as an (FE dimension • K [-dimensional object and building the 

approximation space accordingly (cf. [31, Section 2.3.5]), results in poor accuracy. 
This is due to the fact that the possible savings from reduced problems of smaller size 
seem marginal if achievable at all. In addition high-frequency information can and 
should be exploited for the approximation of low-frequency information, an effect that 
is expected to become more and more important with increasing number of desired 
eigenvalues. In [45], an elastic buckling problem is studied. While the model reduc¬ 
tion is carried out solely/primarily for a linear problem, the eigenvalue problem appears 
only in a second step. 

Here, a new RB space is built from the eigenfunctions associated with the smallest 
eigenvalues at the previously identified parameters. A non-rigorous a posteriori bound 
is then computed by comparison with a reduced space approximation of double size; 
cf. also [43]. Furthermore a component based RB method is studied for eigenvalue 
problems in [42]. 

Very recently an RB method for the approximation of single eigenvalues in the 
context of parameterized elliptic eigenvalue problems has been investigated in [12]. 
The authors derive a bound for the error in the first eigenvalue which is assumed to be 
single. 

The aim of this paper is to develop a model reduction framework for elliptic // 
EVPs. The application scenario we target is the vibro-acoustics of cross-laminated 
timber structures. Here, a parameter-dependent eigenvalue problem in linear elasticity, 
where the input parameters are the material properties of different structural compo¬ 
nents, is to be solved many times during a design/optimization phase. Since the main 
part of a vibro-acoustical analysis is the modal analysis, which not only takes the first 
eigenvalue into account, but all eigenvalues under a certain frequency depending on 
the problem under investigation, the outputs of interest are the K smallest eigenvalues 
with corresponding eigenfunctions. A characteristic feature of the considered /iE VPs 
is the appearance of multiple eigenvalues. In particular, the multiplicities depend on 
the parameters. 

A particular challenge of the considered /iEVPs is the rather large number of out¬ 
puts of interest K , which in our exemplary case ranges from two to twenty. We are 
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interested in approximating these smallest eigenvalues “simultaneously” in the sense 
that a single reduced space is constructed for the variational approximation of the eigen¬ 
value problem and that the individual a posteriori error estimators for the eigenvalues 
use the online components provided offline. This allows us to generate an efficient and 
accurate simultaneous reduced basis approximation. The large number of outputs of 
interest K justifies an increased computational effort by an increased dimension N as 
compared to the standard single output case. In particular, any decent, i. e., sufficiently 
accurate, approximation needs N K. This approach constitutes a significant differ¬ 
ence to the one taken in [12] as our goal is a reduced basis approximation not only of 
one eigenvalue, but of a series of eigenvalues, including eigenvalues with multiplicity 
greater than one. The parameter-dependence of the multiplicity of the eigenvalues con¬ 
stitutes a major challenge in this context and is included into both our analysis and our 
algorithms. 

Furthermore our experiments show that, in a greedy algorithm, it is usually not 
optimal to include the first K eigenfunctions for a particular parameter, neither is it 
advisable to choose the same number of eigenfunctions for different eigenvalues. This 
may be attributed to the fact that the smoothness of the input-output relation can vary 
strongly with the different outputs of interest, i. e., the eigenvalues. We rather sug¬ 
gest to choose maximizing parameters for K different error estimators, as described 
in Sect. 4.1. The reduced approximation should be of comparable quality for a broad 
range of frequencies, although in structural acoustics the accuracy requirements might 
decrease with increasing frequency. Note that, for the application scenario at hand, the 
number of desired eigenpairs is typically in the order of ten for simple components and 
even larger for geometrically more complex structures. 

The main contributions of this paper are the analysis of an asymptotically reliable 
error estimator including the case of multiple eigenvalues and a series of algorithmic 
advancements. Our numerical results demonstrate that tailored greedy strategies yield 
very efficient reduced basis spaces for the simultaneous approximation of many eigen¬ 
values for the considered problem class. 

The rest of the paper is structured as follows: In Sect. 2, we describe the problem 
setting and introduce the reduced basis method for /iEVPs. Sect. 3 is devoted to the 
a posteriori error analysis w. r. t. parameter variations. We also discuss how to evaluate 
the derived error estimators computationally. In Sect. 4, several greedy algorithms are 
presented. We demonstrate the effectivity of our algorithms by numerical examples 
with the application to linear elasticity in Sect. 5. 

2 Problem setting 

2.1 Parameterized eigenvalue problems in computational mechan¬ 
ics 

Let the computational domain Cl C with d = 2,3, be bounded and polygonal. As 
an elliptic eigenvalue model problem, we consider the linear elasticity case. But all our 
results also hold true for more general elliptic systems. Then, the eigenvalue problem 
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in linear elasticity is given by 


—diver = Xpu in Q. 


( 1 ) 


with boundary conditions prescribed as Dirichlet conditions on a closed non-trivial 
subset r D of dQ. and homogeneous Neumann conditions on r)£2 \r D . In addition, the 
linearized stress and strain tensors are defined as 

1 T 

off) = C(p)e(w) and e(u) = -(Vu + Vu ), 

respectively. We set the density p to 1 for simplicity. Furthermore, the set of admissible 
parameters is denoted by 2? C and Li £ /X 5 stands for a vector of parameters. Then, 
C(jtt) denotes the parameter-dependent Hooke’s tensor, which we assume to be uni¬ 
formly positive definite. To this end, let £2 be decomposed into non-overlapping sub- 
domains such that £1 = [J, £2 S . We assume that the material parameters are piecewise 
constant w. r. t. this decomposition. In the isotropic case, the parameters may be chosen 
as Young’s modulus E and Poisson’s ratio v such that P equals two times the number 
of structural components (i. e., subdomains). More precisely, we set p,2s-i = E |n, and 
P 2 .s = v|n. in this case. The anisotropic case is treated analogously. 

Let the bilinear forms a(-,-;ju): (H l (Q.)) d x (//' (£2)) d — y R and m(-. •): (L 2 (£2)) d x 
(L 2 (£2)) d —>• K be given by 

(u,v) a(u,v;/i) := / C(ju)e(i<): e(v)dr 

Ja 

and 

(u,v)^m(u,v) := (u,v) L 2 (a) := J^u-vdx. 

Note that a(-, -;p) depends on the parameter vector u whereas m(-, •) and £2 do not. 

Remark 2.1. The equations of linear elasticity are used as a model problem as we 
are interested in the applications of vibro-acoustics. However this does not pose any 
restriction to the theoretical results shown in the following. Thus we could replace 
a(-,-,p) by any elliptic bilinear form. 

Let V C {v £ ( H l (£2)) d \ vf r) = 0} be a fixed conforming finite element space of 
dimension jf. Then, the discrete variational formulation of (1) reads as: Find the 
eigenvalues A(ju) £ R and the eigenfunctions u(ju) £ V such that 

a(n(ju),v;ju) = A(p)m(«(p),v) VvSV (2) 


for given ju £ 2P. We assume that the eigenvalues are positive and numbered as 

0 < Ai(ju) < ... < Xjy(p). 

The corresponding eigenfunctions are denoted by «,-fp) £ V for i = 1,.... ,/L with the 
normalization 

m(ui(p),Uj(p)) = Sjj for 1 < i,j < jV. 
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In the present context, the error of the finite element solution is assumed to be very 
small. This is achieved by a fine mesh size leading to a large dimension V. The 
discretization error analysis can be found, e. g., in [1-3]. 

Let L > 1 be the number of distinct eigenvalues of (2). For multiple eigenvalues, 
we use the standard notation from [3] and denote the lowest index of the /-th distinct 
eigenvalue by kj and its multiplicity by qt, i = 1,... ,L. We write := {£,, + 

qi — 1}. (Here and in the following, the dependency of the index notations on /i is 
suppressed as it is always clear from the context.) The corresponding eigenspaces are 
denoted by 

Ut(H ): = span{ u kj (ju),..., u k . +q ._ x (ju)} . 

Now, the goal is to find a computationally inexpensive but accurate surrogate model 
that can be used in the many-query or real-time context. 

2.2 Model reduction 

We consider a variational approximation of the /xEVP in an Af-dimensional reduced 
space 

V N := span{£„ | n= CV, (3) 

N <C ,/T. As a matter of fact, the choice of Viv highly depends on the algorithmic 
methodology. Several (snapshot-based) possibilities are investigated in Sect. 4. 

Now, the “reduced eigenvalue problem” reads as 

(u re d(At), Aj. e d(ft)) € V/v x IR, tt(w re( i(lt), v,) = A rec i(ft)/rz(M re d(lt), v) V v € Viv (4) 

for given ji £ £?. Let us emphasize that all eigenpairs of interest are approximated 
in the same space Vjv- As before, we assume a numbering dred,;(M)i i = 1,... ,N of 
the “reduced eigenvalues”. The minimum-maximum principles guarantee that Xi(fi) < 
Ared.Kft) f° r i = 1,...,1V; see [3, Sect. 8]. Note that the multiplicity of the finite el¬ 
ement eigenvalues is not necessarily reflected in the reduced basis eigenvalues. The 
corresponding eigenfunctions are denoted by M re d.,(/i) £ Vi v for i= 1,..., N, again with 
the normalization 


m (n re d,;' (M)) “red, f (M)) — for 1 < i. j £ N. 

In practice, as mentioned before, one is only interested in the first K eigenvalues for 
any chosen parameter. We expect that the dimension N required to achieve a certain 
accuracy will depend not only on the smoothness of the parameter-dependency of the 
/xPDE but also on the number of outputs K. 

In the present context, a(-, ;/r) is affine w. r. t. the parameter li, i. e., 

Q 

a(u,v\n) = Yj (5) 

9=1 

for suitable parameter-independent bilinear forms a q : (II 1 (Q)) d x (II 1 (Gi)) d R and 
coefficients © 9 : £? —>■ R, which are readily derived from the constitutive equations. 
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For instance, we have two terms per subdomain in the isotropic case. This leads to a 
fast online evaluation as the cost of the assembly of the parameter-dependent reduced 
systems (i. e., matrices in M. NxN associated with (4)) is independent of ,jY . Note that 
the expansion (5) will also be exploited for an online-offline decomposition of the error 
estimators. 


2.3 Model reduction by proper orthogonal decomposition 

Before turning to the development of greedy methods and a posteriori error estimators 
for /iPDEs, we illustrate the potential of model reduction techniques in the context of 
parameter dependent eigenvalue problems with multiple output values. 

A common technique for model reduction is the proper orthogonal decomposition 
(POD) [21, 33, 34], which yields the best possible reduced space (in the sense that, for 
a given series of snapshots, the projection error w. r. t. the Zr-norm is minimized). To 
this end, let S C V be a set of snapshots generated by solving the /iEVP (each time for 
K eigenfunctions) for all parameters in a sufficiently large training set S[^®. Then, in 
the definition (3) of the reduced space V/v, orthonormalized functions {Ci, - - - - Cn} C 
spanl-S} are selected such that 

Lllv-rMI^) 

vGS 

is minimal. Here, ITv is the /.--orthogonal projection to span {£] n}- This essen¬ 

tially amounts to assembling the #S x #5-correlation matrix of the snapshots in S w. r. t. 
(',')l 2 (£ 2 ) anc * finding its N largest eigenvalues and corresponding eigenvectors. For 
a detailed description of the usage of POD methods in the present context, see, e. g., 
[21, 33, 34], 



Figure 1: Convergence of POD methods (0 < N < 200) for the described /iEVP with 
different numbers of outputs of interest: Average relative errors in the eigenvalues 
A|,.... Xk for K = 4 (left) and K = 7 (right) 

Fig. 1 illustrates the convergence of the POD method for K = 4 and K = 7. (The 
details of the underlying numerical experiment are elaborated in Sect. 5.) On the one 
hand, the results show that it is in principle possible to construct one single reduced 
space that effectively captures the parameter-dependent behavior of the first K eigen¬ 
functions simultaneously. On the other hand, it is evident that the RB dimension N 
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required for a certain accuracy increases with K. More precisely, the asymptotic decay 
of the error is approximately C 4 e~ 0 • 0513 ■ A, for K = 4 and C 7 e _004777V for K = 7 for 
some constants C 4 and C 7 . From the cost point of view the POD is quite expensive, 
and thus we focus on computationally efficient greedy strategies in combination with a 
posteriori error bounds. 


3 A posteriori error estimation 

In this section, we will establish a posteriori error estimators for our output quantities, 
in this case the eigenvalues. It is important to determine such estimators in order to 
find out which basis functions should be selected by the greedy method. In particular, 
their computational evaluation must only depend on the basis size N but not on the 
dimension of the finite element space JV. 

To this end, we first derive error bounds that still depend on the finite element 
eigenvalues, and in particular on their multiplicities. Then, Sect. 3.2 is devoted to a 
computable approximation yielding the desired error estimators. 

In Sect. 3.3, we recall a standard online-offline decomposition. 


3.1 Error bounds 

Let the parameter-dependent energy norm be defined as |M|^ y : = a(■, ■; ,U ) 2 . 

In addition to a parameter dependent norm we are using a parameter indepenent 
norm defined as ||-HjQ-v :== a(v) z : = a(-, •-A) 2 • For a linear functional r : V —> R, we 
define the corresponding dual norms by 


171 ;t" 


, := sup 


r(v 


and 


O^vGV \\ v \\n;V 


r(v) 


(i-y 


respectively. 

The analysis and the practical implementation employ different error representa¬ 
tions, namely the so-called reconstructed errors w. r. t. the bilinear forms a(-, •; jJL) and 
a(v). Using the residual 


v !-)• n (v;ju) : = a(u re d,;(£ 0 ,v;/f) - -W,• (ju) m( M red,;(jU), v) 
for i = 1,... ,N , we define e,-(ju) G V and e, (jtl) € V by 


a(e,-(ju),v;ju) =n(v;ju) Vv€V 


(6) 


and 

a(ei(n),v) = ri(v;n) Vv€V, (7) 

respectively. In particular, ||l(-;m)II m; v' = Ik/(M)ll /1; v and \\ r i(-^)\\fi-v' = \\ei(^)\\fi-y 
For any Ll £ assume that g(!-L) > 0 is a generalized coercivity constant such 
thatg(jtl)a(v, v) < a(v, v\n) for all v € V. (Technically speaking, g(ju) is the parameter- 
dependent coercivity constant of a(-, -;/f) w. r. t. ||-||^.y.) This implies 

ll r ll fl;V' — S(^) ~ IMI/t;V' 18) 
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for any r £ V' and 


(9) 


l|v|| A; y<^r 2 l|v|| M; y 

for any v € V. 

We are now ready to prove the error bounds. The following theorem, combined 
with the computational/algorithmic aspects in Sect. 3.2 and Sect. 4, generalize the re¬ 
sults of [24, 30-32] for the case of multiple eigenvalues. 


Theorem 3.1. Let 1 < i <L such that k + g; —1 <N. For j = \,... ,q t , set 


d-kj+j— l (m) • 


min 

jY>l>kj+qj— 1 


A/(AO ^red.ki+j— 1 (M ) 

MaO 


( 10 ) 


Then, 


0 ^ Ared, ki+j— 1 (M ) A^. (M ) — 


jujv/ ^ ^i+y'-d'iA 4 ) j, ; y/ 

<4,-+;-l (AO dki+j- l{H) 2 Vhi 


(ID 


Proof. Fix /i G , 1 < i < L and 1 < j < q\. Let M re d,Ar,+y— l (A 1 ) = L/=i «/«/(At) and 
1 (M) = I/=i j3/«/(At)- B y (6). we find 


ft 


A/(At) A rec [, t',+7 — 1 (At) 

4(4 


Therefore, we get 


^ i+ /-i(-;At) 


2 

/r;V' 


> 

> 


2 (hed,kj+j-l(P)\ 
a l 


I—l 


M At) 


A 


A/(At) 


E «/ 


2 /A/ (M ) ^red, &,+ 7 — 1 (M ) 


A, (ai) 


A/(At) 


4+2-1 (At) 2 E «/%(/t)- 

l^>kj-\~C[i — 1 


( 12 ) 


Using the fact that af = 1, A/(/x) < A^. (/i) for l < hi +qi — 1 , we find for the 
difference between approximated and detailed eigenvalue 


• — Aj- ec i_ 7 — 1 (fi ) A&. (/i) 


^(Wred,&;+7— 1 (M ) ? M red, ki+j— 1 (M ) » M ) A&- (/X ) 

Ea/ 2 A/(At)-A A .,(At) 

;=t 

E a, 2 (A, 70-At, 70)+ E «/ 2 (4/(At)-4(At)) 

K*i+? — l l>ki+qi~t 

E 4 (A/ (At) — Aa 7 (At)) ■ 

l>k t +qi -1 
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< 


(13) 



From this we obtain two upper bounds for AA^. The first one follows trivially from the 
fact that A^.Qll) > 0 and (12) 


AA*,. < E a, 2 A i{n) < 

l>kj+qi -1 


r ki+j~ l('iM) 
d/q+j-l (fO“ 


(14) 


The second bound is based on the Cauchy-Schwarz inequality and on Young’s inequal¬ 
ity. In terms of, 

(A/(jU) — Aa- ; (/X)) 2 < (1 + e)(A;(/l) — Ared^+j-ltrt)) 2 + (1 + -)(AAfc,.) 2 

for e > 0, we get from (12) and (13) 


AAt, = 


< 


E y/k{v)aiy/Mn) 




< 


< 


< 


V >£;+?;-1 

1 


E «/%(/*) 

i / 


dt 


+j- 


\\ r k i +j-l(''f i )\\n-V'^j( l+£ ') ll r ^+f-l(-;M)||^.y/ + ( 1 + g) A A^ ^E 


1 


4,+y-l(M) 

1 , 
dki+j— l (a 0 

lk-+7-i(*^ 


k+j-i(-;/x)|| M;V 


/ T 1 AA, 2 

I (!+ e ) || r *i+7-i(-;M)||“ V '/ + (l + -)■ 


e 4 ,(a0 


2 

tiT' 


\ 


1 K+j-i (•;/*) 
l + e + (1 + -) 


2 


® dkj+j-\{[l)^^kj{P-) 


Setting e = y ' IJ . .., .A 1 ','' gives the upper bound in (11). The lower bound follows 
directly from [3, Sect. 8], 

□ 

Besides the generalization to multiple eigenvalues, let us point out that our bounds 
are sharper than the ones, e. g., in [32, Prop. 1], as the lowest order term in (11) is of 

the form rather than 


%+j -1 


(ju) 


dkj+j-lW 2 


Note that the error bounds in 


Thm. 3.1 still depend on the finite element solution via the eigenvalues A/(jU) in (10). 

Remark 3.2. It is also possible to give an upper bound for the eigenvectors by replac¬ 
ing d^+j- 1 (ft) by dki+j-1 (aO defined as 


dk+j-\(u) := min 

<jY'>l>ki+qi— 1 V l<ki 


^/(AO ^red.ki+j— 1 (AO 

A/(jU) 
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Using now IT/: V —> Ufip) as the orthogonal projection w. r. t. the Lr-inner product, we 
define v := n,(M re d,* i+ y-i (fi)) = L^>/> 4 ,-+ 9 ,-i v/<i, «/ M /(M) and give the upper bound 
as 


^red.A, — / ] (/-0 



2 


£ aw (ju) 

<S^>l>ki+qi— 1 V l<ki 


£ a, 2 A/(Ai) < 

■jV>l>k ! \ q, I V /<£, 


r ki+j- l(‘>M) 
dkj+j— i (aO“ 


3.2 Error estimators 


We now derive approximate error bounds that are computable in the sense that they do 
not depend on the finite element solution. To achieve this, it remains to approximate 
di, which may be interpreted as a measure for the relative distance between neighbor¬ 
ing/adjacent eigenvalues, particularly to decide which of the indices to exclude from 
the minimum. We point out that the dimension of the (detailed) eigenspace is not ac¬ 
cessible. The application scenario we have in mind features multiple eigenvalues with 
their multiplicities depending on the parameter. It is therefore impossible to determine 
the structure of the spectrum (i. e., the indices k, or the index sets ./fi) a priori. 

Recall that the first K eigenvalues are the output quantities of interest. Assume that 
the reduced basis method converges in the following sense: For p £ & and 1 <i<K, 


Ared.KfO A,(jli) for N —> ,A . 


In particular, X K( i j(p ) —»• X^fip) for N — > JV for j € Jjj. 

Given the eigenvalues A re d, ( (jtl)> i = 1 of (4), we replace A/(jtt) in (10) by 

Ared./(jU.) and approximate by 


^red a ■= { 1 < j <K + x: 


Ared, j ( P ) A ret | i (p ) 
Ared, /(M) 


<£2, 


for a chosen “tolerance” > 0 and with t as the difference between the index of the 
first eigenvalue after the multiplicity of the Zf-th eigenvalue and the Zf-th eigenvalue 
itself. In the case that we know a priori the maximal multiplicity of all relevant eigen¬ 
values for all parameters, we set r equal to this value. Otherwise we select it adaptively 
during the initialization phase of the greedy method. More precisely, we start with r = 1 
and increase it by one as long as K +1 £ d, /• Thus #J^ e d. / will be our best guess for 
the multiplicity of the eigenvalue to which A re d.if,U) converges. Then for 1 < i < K, 


di{p) := 


min 

/^J^red, / 

K+v>l>i 


Ared, lifi-) Ared ,i{P) 
Ared. / ( M ) 


(15) 


is the relative distance of A ret i , (/d) to the reduced eigenvalues that are further away than 
the chosen tolerance e^. The adaptive selection of r guarantees that even for i = K and 
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multiple eigenvalues dj(fl) is easily computable and does not severely underestimate 

di(n). 

Finally, since we are looking for an asymptotic estimator for the relative error in 
the eigenvalues which is cheaply computable in the online-phase, we neglect the higher 
order term in ('ll). In addition, the parameter-dependent norm H'll^.y, is replaced by 
the parameter-independent norm ||-||^,y/ by (8), which introduces an additional factor 
g(p)^ 1 . To summarize we can state the following corollary: 

Corollary 3.3. Let i = 1 ..... AT and A re d,<(M) —> Xfip) for N — > jV . Furthermore let 
JLied.i be defined as above and the distance between neighboring eigenvalues dfi/J,) be 
given as in (15). Then the error estimator given by 


H ^ rji(n) := 


g(ju)-rf/(ju)-A re d,i(M)' 


( 16 ) 


is asymptotically reliable in the sense that 


0 < Wj-tOO A,,(At) < 
hi 

with C tending to one as N tends to Jfi. 

Note that the approximation in (15) is, in general, less accurate for i = K. This 
is because the space Vn is built to approximate well the K outputs, but for the Af-th 
estimator we need the (AT + t)-th outputs with r > 1, which are approximated only 
roughly. The tolerance Ex has to be selected such that it reflects the desired accuracy 
of the RB approximation. 

3.3 Online-offline decomposition 

All error estimator contributions may be decomposed as already outlined in [24], Let 
(Cn)i<n<N be the orthonormal basis (w. r. t. m(-, •)) of Vn- For 0 < q,p < Q let A chP £ 
M. NxN with Anfm := a(^n,^m) for 1 < n,m < N where 

a&, v )= a q(£n,v), VveV, l<n<N, l<q<Q, (17) 

and 

a(4,?,v) = m(£„,v), VveV, 1 <n<N. (18) 

In the following, we identify the function u m ij(p) £ Vn and its vector representation 
w. r. t. the basis ( Cn)i<n<N such that (wred.KA 1 )),, denotes the n-th coefficient. Then, 
given a reduced eigenpair (u Te d,, (jtl). Ar e d,(jtt)), we have the error representation 

Mm) = E E 0 ?(M)(Mred,i(M)) n 4« 9 -^red,i(M) E (“red,4n 

n= 1 q =1 n =1 
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by (7). Consequently, the main contribution of decomposes into 

lki(-;/i)llA;V' = L E E E Kd > i(M))„(Mred,i(M)) m ©?(M)©p(M)^ 

h= 1 m=l q= 1 /?=1 

N N 

+ ^-rcdE E ( M red,i(M))„ ( M red 

n= 1m=l 
N N Q 

— 2A re d, i '(/j) y E E ( M red ,i(n)) n ( u red.i(H))„, At.m- 

«=1 m=l ty=l 

We recall that only a single reduced space is built for the approximation of all 
eigenvectors simultaneously. Thus the above decomposition uses the same offline in¬ 
gredients for all 1 < i < K. In particular, the number K of desired eigenpairs does not 
directly influence the complexity (only via the reduced space dimension N). 

4 Algorithms / Basis construction 

In this section, we present different greedy strategies that employ the error estimators 
of Sect. 3 to build the reduced space in (3). The advantage as compared to the POD 
method motivated in Sect. 2.3 is that only relatively few finite element solutions of the 
jUEVP need to be computed. 

Since we use a single space for the approximation of multiple outputs, we have 
several natural possibilities which are investigated in Sect. 4.1. Sect. 4.2 is devoted to 
an extension that takes into account multiple eigenvalues. In Sect. 4.3, a remedy for 
the potential unreliability of the error estimators for small N is discussed. 

4.1 Greedy selection of snapshots for single eigenvalues 

Recall that the K smallest eigenvalues are the quantities of interest, where K is typically 
2 — 20 for our application scenario. In principle, given a reduced space, one could try to 
identify a suitable Ll £ '¥ and then include the first K eigenfunctions for this parameter 
value. (In each greedy step, this would require the detailed FE solution of (2) for 
one parameter only.) However, numerical studies clearly show that this naive choice 
is far from optimal as the generated reduced spaces tend to be much too large. This 
is because the errors in the individual eigenvalues and eigenfunctions are only very 
weakly correlated, if at all. There are at least the following two much more natural 
options. 

Let a sufficiently rich training set Strain C £? be given. Then, in Alg. 1, the indi¬ 
vidual argmax for each 1 < i < K is chosen separately. In contrast, Alg. 2 chooses 
only one single argmax. Note that both Alg. 1 (line 7) and Alg. 2 (line 6) require the 
evaluation of all error estimators at all parameters in E tra i n to determine the choice of 
fi. This does not lead to large computations since the calculations are only performed 
with the reduced space of size N , such that we obtain K reduced eigenpairs for any 
u £ S tra i n , see a ls° Sect. 3.3. However, Alg. 1 (line 10) and Alg. 2 (line 8) require also 
finite element solutions which then determine the reduced basis space. 
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Algorithm 1 Multi-choice greedy 
l: for i = 1,. .. ,K do 

2: C 

3: end for 

4: N-h- K 

5: while N < AW do 
6: for i = 1,..., K do 

Mmax, < ar g rnax ( ues traln 
8: if max,/) £tol then 

9: A^/V+l 

10: £/v f- Ui(jX mm i) (orthonormalized) 

li: end if 

12 : end for 

13: if max MeHtrain) i</<jf t]> (ju) < e to i then 

14: break 

15: end if 

16: end while 


The multi-choice variant rests on the intuition that the individual eigenfunctions 
can/should be approximated separately. In contrast, the single-choice variant takes into 
account that the approximation power of eigenfunctions to large eigenvalues can be 
exploited also for eigenfunctions to smaller eigenvalues. 

During the greedy procedure, we orthonormalize the selected basis functions. Not 
only does this yield small condition numbers of the reduced systems; it is also benefi¬ 
cial for the special treatment of multiple eigenvalues described in the next section. 

In Alg. 1 (line 10) and Alg. 2 (line 8), an orthonormalization is performed. For this 
purpose, let ITy : V —► Vjy be the /^-orthogonal projection to the current reduced space. 
For a snapshot candidate £ € V (i. e., one of the eigenfunctions chosen as described 
above), we compute £ := £ — FIyvC- Then, if ||£|| is sufficiently large (> £ pro j), the 

“new contribution" ? is included in the reduced basis; see also Sect. 4.2. 

4.2 Extended selection for multiple eigenvalues 

In case of multiple eigenvalues, the greedy method needs to be modified as follows. 
Assume an index 1 < i < K and a parameter p. have been selected by means of the 
eigenvalue-based estimators ( 77 ,)f = 1 ,in Alg. 2 (line 6), or several parameters fi 
have been selected in Alg. 1 (line 7), such that the span of uj(fl) is to be included in the 
reduced space. 

However, a large value of ) merely indicates that the corresponding (fine) 
eigenspace Uj(jl) contains functions that are badly approximated by the current re¬ 
duced space. Nevertheless the eigenspace might also contain other functions that are 
already well approximated. Consequently, if the detailed eigenvalue associated with 
a chosen snapshot has multiplicity greater than one, we aim to add all the eigenfunc¬ 
tions for the multiple eigenvalue, except the ones which are already approximated well 
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Algorithm 2 Single-choice greedy 
l: for i = 1,... ,K do 

2: C 

3: end for 

4: N-h- K 

5: while A < A max do 

6: (/4 max i (max ) <- argmax^gg^ 

7: A -S— A+ 1 

8: CAT^-w/maxCMmax) (orthonormalized) 

9: if max M6Straini1 <K/f ?7/(/r) < e to i then 

10 : break 

11: end if 

12 : end while 


enough. A motivation for exploring the whole eigenspace for multiple eigenvalues is 
to guarantee that we take the correct eigenvalue/eigenfunction, since we cannot ensure 
that the indexed eigenvalue/eigenfunction in the reduced space is the same as in the 
detailed calculation. This is due to the fact that there is no prescribed ordering for the 
eigenfunctions corresponding to a multiple eigenvalue. 

As for the definition of cl, (fl) one has to compute a sufficient number K' > K of 
eigenfunctions of the finite element problem fiEVP (2) such that X K i (fl ) / Xk (fi) > 
1 +£\. Then, lines 9-10 in Alg. 1 are replaced by: 

for all 7 > 1 with Xj (.Umax.:) Xj (.Umax./) |/ Xj (/-/max./) T do 

if IT/(/-tmax,/) Fly: llj (ftmax,/) IIr 2 (tt) — ^proj then 

A^-A+l 

Cn ■<— M/^/tmax,,) (orthonormalized) 

end if 

end for _ 

Analogously, lines 7-8 in Alg. 2 now read as: 

tor all j > 1 with \Xj(^i max ) — Xi max (n max )\/Xi miix (n 

max) ^ e A do 

if !T/(M m ax) n NUj (ji rmx ) |/_ 2 |Qj y £proj then 
A A+ 1 

Cn ^ Uj(n m ax) (orthonormalized) 

end if 

end for _ 

Here, Fly : V —> Ky denotes the /..--orthogonal projection. The parameter e pro j is a 
small tolerance that prevents the selection of functions that are already approximated 
sufficiently well. 

4.3 Initialization of the greedy method 

In our calculations, we need an error estimator for the AT-th eigenvalue. For the com¬ 
putation of this estimator, we need a rough approximation of the (K -|- t)-th eigenvalue. 
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In order to ensure that our reduced space has the ability to roughly approximate this 
(.K + t)-th eigenvalue, we use an initial approximation space in which we include the 
corresponding components. We suggest to include components using the proper or¬ 
thogonal decomposition method described in Sect. 2.3 (with N = /V m it) applied to a 
small number of snapshots. Here, the snapshots S C V are associated with a training 
set typically of size 2 P , taking into account the extension described in Sect. 4.2. 
This initial approximation space of dimension which is constructed as an initial¬ 
ization step for the greedy algorithm, should be sufficiently large as the reliability of 
the error estimators analyzed in Sect. 3 can depend on the dimension of the reduced 
space. To make sure that we are able to calculate and to approximate the (K + r)-th 
eigenvalues, we chose our A^ n i t to be at least (K + r) times a factor > 1.5. 


5 Numerical results 

In this section, the performance of the proposed algorithms is illustrated by numerical 
examples, in two and in three dimensions. For the two-dimensional calculations we 
use plane strain elasticity while for the three-dimensional simulations we use linear 
elasticity. The implementation is performed in MATLAB based on the RBmatlab [9] 
library. We investigate the individual components and highlight their benefits in several 
steps. 


5.1 Preliminaries 


First, in Sect. 5.2 to Sect. 5.5, we choose £2 as rectangle of size 3.0 x 1.0 with Dirichlet 
boundary on the left and on the right. Let £1 be split into three subdomains of size 
1.0x1.0. The material parameters E and v used for these subdomains are in the range 
of 10 — 100 and 0.1— 0.4, respectively; we have P = 6 and Q = 6. We choose a uniform 
random sample of size 10,000 as set of training parameters E tra j n C . To evaluate the 
errors, another sufficiently rich set of parameters E test C is used of size 1000. For 
our initial space we choose A; n it < 40, depending on the desired number of eigenvalues 
K. We always report the average errors of the reduced approximations given by 


1 


#“test 


£ 

M^—test 


^red ,i(M) 

MaO 


and comment on the standard deviation at the end of Sect. 5.4 in Remark 5.1. 

For the generalized coercivity estimate, we exploit the affine decomposition of the 
bilinear form and set 


g(M) := 


. ®„Gu) 

min • 

«=!,-,0 ©, y ( A 0 


(19) 


We emphasize that g(fi) merely relates the bilinear forms and a(-,-); a coer¬ 

civity estimate for a(-, ■ ; ,Ll) itself is not required in the present context. Note that (19) 
indeed yields an admissible parameter-dependent constant provided the bilinear forms 
a q {•,•) in (5) are positive semi-definite and the coefficient functions © f/ (-) in (5) are 
positive; see, e. g., [29, Sect. 4.2.2]. This is true for our application. Better results (i. e.. 
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a larger lower bound) could be obtained by the more expensive successive constraint 
method [19]. In the present setting, the estimate (19) is typically smaller than the exact 
solution of the corresponding generalized eigenvalue problem by a factor ranging from 
0.7 to 0.98. 

5.2 Extended selection vs. non-extended selection 

We first illustrate the necessity of the extended selection for multiple eigenvalues. 
Fig. 2 shows the behavior of a POD method with (left) and without (right) the extended 
selection described in Sect. 4.2 for the first two eigenvalues (K = 2). Fig. 3 shows 
the same comparison for the greedy method (Alg. 2). For both the POD method and 
the greedy method, we observe that in the variants without extension the convergence 
for the second eigenvalue becomes slower after a certain number of basis functions 
has been included. In contrast, the extended selection yields convergence curves that 
approximately coincide. 



Figure 2: Error decay for the eigenvalues with the POD method: extended (left) vs. 
non-extended (right). First row: training size 10,000; second row: training size 1,000. 

The shortcomings of the non-extended methods may be explained by the fact that 
the second eigenvalue has multiplicity two for certain parameters and in these cases, 
for the multiple eigenvalue, the correct eigenfunction is not necessarily chosen. Note 
that the effect is more significant for a smaller POD training size (Fig. 2, second row) 
as it is less likely that all directions of an eigenspace are present in the snapshot set. 
The convergence of the second reduced eigenvalue possibly improves drastically if, 
incidentally, the missing component is added during the greedy method. 
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Figure 3: Error decay for eigenvalues with the greedy method: extended (left) vs. non- 
extended (right); training size 1,000. 


5.3 Multi-choice vs. single-choice greedy method 

Here, we illustrate the benefit of Alg. 2 in comparison to Alg. 1. In Fig. 4, for K = 4, 
one can see that with Alg. 1 (left) the convergence behavior varies over the course 
of the greedy method while with Alg. 2 (right) all desired eigenvalues exhibit similar 
convergence. (This also holds true for the errors in the eigenfunctions not shown here.) 



Figure 4: RB error decay: comparison of Alg. 1 (left) and Alg. 2 (right) for K = 4 

The poor convergence of the third eigenvalue only improves rapidly at N ~ 170, 
after the other three eigenvalues have reached an accuracy in the order of the target tol¬ 
erance, and thus the algorithm only chooses EV 3. This effect (namely an imbalanced 
resolution of the relevant eigenspaces during the greedy method) is directly related to 
the inappropriate a priori assumption of Alg. 1 that roughly the same number of snap¬ 
shots corresponding to the first K eigenvalues should be included in the reduced space. 
At this point it should also be noted that in general Alg. 1 creates a larger RB space 
than Alg. 2 as soon as more eigenvalues have a poor convergence. 

To further illustrate the behavior of the single-choice greedy method, in Fig. 5, we 
report the accumulated numbers of chosen eigenfunctions corresponding to At,..., A/f 
over the course of Alg. 2 for K = 4 (left) and K = 7 (center), as selected by the error 
estimators in line 6. The reason for the greedy algorithm not selecting any eigenfunc¬ 
tions before a basis size of 40 is that this is the size of our initial space. The respective 
error decay for K = 7 is depicted in Fig. 5 (right). Note that the good convergence 
(in particular, similar rates for all outputs of interest simultaneously) is achieved by 
a rather uneven distribution. The diagrams indicate that, for both values of K, larger 
eigenvalues as well as possibly double eigenvalues are preferred by the algorithm. This 
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Figure 5: Left and center: accumulated numbers of chosen eigenfunctions over the 
course of Alg. 2 for K = 4 and K = 7. Right: error decay for K = 7 


and the fact that, although fewer eigenfunctions are included for the smaller eigenval¬ 
ues than for the larger ones, but nevertheless the error decay is equal, mean that the 
eigenfunctions corresponding to larger eigenvalues are effectively used to approximate 
the ones corresponding to smaller eigenvalues. 

5.4 Effectivity of the greedy method 

In this section, we investigate the performance of the greedy method in more detail. 
For this purpose, we also consider the effectivity numbers Ji, 1 < i < K, of the error 
estimators and its maximal ratio R defined by 

1 y r]i(n)-ktii) R ._ max,- |. K y, 

# £test ^ est Ared,;(jU) - A/(ju) ’ ' min, = i K J, ' 

As already mentioned, the estimators derived in Sect. 3 are of asymptotic character and 
therefore generally not reliable for small N. To prevent a misleading selection of basis 
functions in the first few greedy steps, the initialization described in Sect. 4.3 is used 
to generate an initial basis. 

Fig. 6 shows the error decay (left), the effectivity numbers of the a posteriori esti¬ 
mators (center) and the accumulated index counts (right) for K = 5 with and without the 
initialization. In this case, a similar convergence is achieved for both algorithms, and 
the index count plots also shows a similar behavior. In the preasymptotic range, we 
observe a difference in the effectivity numbers. Without initialization these numbers 
possibly depend sensitively on the selected snapshots. While this does not influence the 
overall performance for K = 5, for K = 7 we do get extremely poor results if we start 
directly with the greedy algorithm. This is caused by the fact that the approximation 
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Figure 6: Comparison of Alg. 2 without (top) and with (bottom) the initialization de¬ 
scribed in Sect. 4.3 for K = 5 for a selected case in which the greedy algorithm without 
initialization does not fail. RB error decay (left), effectivity numbers (center) and ac¬ 
cumulated index counts (right) 


of di(fi) by dj(fi) is then not reliable. Thus we always include the initialization step in 
our adaptive algorithms. 




Figure 7: Effectivity numbers of the estimators for K = 4 (left) and K = 7 (right) 

In our experiments, the described initialization always prevents the effectivity num¬ 
bers from having jumps and leads to good convergence of the greedy methods. For 
instance. Fig. 7 shows the effectivity measures corresponding to the error curves from 
Fig. 4 (right) and Fig. 5 (right). The effectivities are virtually constant and close to¬ 
gether which is reflected in a small value of R. This is of crucial importance for the 
performance of our Alg. 2. In all our settings R is below five, e. g., R = 3.41 for K = 4. 
Note that for K = 4 and K = l, the same eigenvalues show similar effectivities. A high 
effectivity ratio R possibly leads to an oversampling of the eigenfunctions associated 
with the indices of a high effectivity and thus a loss in the performance. At this point 
although our error estimators are for eigenvalues, we want to show that also the effec¬ 
tivities for the eigenvectors are constant and close together. To do so we depict the 
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results in Fig. 8. 



Figure 8: Effectivity numbers for the eigenfunctions for K = 4 

After having demonstrated the performance of the single components of our algo¬ 
rithm, let us compare the results of our greedy method using the error estimator and the 
best components with the convergence of the POD method; cf. Sect. 2.3. Comparing 
the error plots in Fig. 1 with the ones in Fig. 4 (right) and Fig. 5 (right), we see that 
we achieve very similar convergence behavior. In particular, the error curves of our 
simultaneous reduced basis approximation for the individual eigenvalues are similarly 
close to each other. Moreover, the accuracy reached at N ss 200 differs only by a factor 
of roughly ten. We recall that the POD method uses the full training set (namely 10,000 
finite element solutions in this case which leads to a computation time of over 10 h) to 
reach this accuracy while the greedy method only needs a couple of hundred detailed 
simulations and the evaluation of the estimator which leads in this case to a compuation 
time of 6-7 h. It should be noted that this gap in computation time between POD and 
Greedy increases further with the complexity of the detailed solution. 

Let us emphasize that the bounds from [32], i. e., c/,(u) 2 in the denominator of (16) 
instead of di(fi ), lead to a large ratio of the maximal and minimal effectivity value and 
thus to poorer results in the multiple output case. Highly different effectivity numbers 
result in an over-selection of eigenfunctions associated with the largest effectivity num¬ 
bers and thus in a performance loss, hence, to a much less attractive greedy algorithm. 


Remark 5.1 (Error evaluation). For completeness, Fig. 9 shows a convergence plot 
including the standard deviation for K = 6. In the semilogarithmic plots, one can see 
that the standard deviation is always in the order of the (relative) discretization error 
itself 

5.5 Speed-up 

The speed-up was calculated serially using MATLAB on a Mac laptop; the standard 
routine eigs, which is based on ARPACK [22], was used for solving the eigenvalue 
problems. We used linear finite elements for the discretization space V. With our 
greedy method as introduced above, a significant speed-up in the computation of eigen¬ 
values can be achieved, as is shown for the settings of K = 2 to K = 7 eigenvalues in 
Table 1. Here, the calculation of the detailed solutions takes in the range of 3.5 to 
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Figure 9: RB discretization error for K = 6 with standard deviation (as unidirectional 
error bar) for EV 2 to EV 5 


3.6 seconds, while the calculation of the reduced solution is possible in 0.021 to 0.078 
seconds, resulting in a speed-up of 140 to 43. The higher the value of N, the longer the 
reconstruction time, but in this case the increase is approximately linear in N. 

Moreover, it should be noted that the more accurately the detailed solution is cal¬ 
culated, the more expensive the detailed calculation becomes while the cost for the 
calculation of the reduced basis solution will stay in the same range, such that we 
would achieve even higher speed-ups. In computations of practical relevance, the de¬ 
tailed and the reduced accuracies have to be adjusted as it is described in [44]. Here we 
are mostly interested in the performance of the RB algorithm, and thus we work with a 
fixed moderate finite element resolution of 15402 DOFs. 


K 

2 

4 

7 

N 

50 

100 

150 

200 

50 

100 

150 

200 

50 

100 

150 

200 

Detailed solution 

3.5 

3.5 

3.6 

Reduced solution 

0.025 

0.028 

- 

- 

0.037 

0.043 

0.050 

0.055 

0.057 

0.067 

0.075 

0.082 

Reconstruction 

0.009 

0.011 

- 

- 

0.009 

0.012 

0.016 

0.019 

0.010 

0.013 

0.018 

0.021 

Speed-up 

140 

125 

- 

- 

94 

81 

70 

63 

63 

53 

48 

43 


Table 1: Timings for the detailed solution and the online calculations (reduced solu¬ 
tion including error estimation; reconstruction of the finite element solution from the 
reduced solution) in seconds and speed-up numbers 

As can be seen in Table 2, the computation times for the error estimators (rj) as well 
as for the required offline components for the error estimators, i. e., solution of (17) and 
(18) (“Assembly”) and computation of A, increase for increasing values of N. In the 
case of A, the increase is approximately linear. Note that these longer computations 
will only have to be performed in the offline phase and will not have any impact on the 
computation times for the online phase. The computation of g(fi) as defined in (19), 
which is necessary for the error estimator, takes 0.0042 seconds. 


21 











































N 


A 

Assembly 

50 

0.0028 

1.3672 

6.076 

100 

0.0037 

2.7640 

6.156 

150 

0.0040 

4.1144 

6.272 

200 

0.0046 

5.5995 

6.312 


Table 2: Timings for single components of the offline phase in seconds 


5.6 Wall-slab configuration 

In this section, we show the ability of the newly developed reduced basis method to 
approximate multiple eigenvalues in a two-dimensional wall-slab configuration with 
a thin elastomer layer in between. The domain shape is an L-shape with three non¬ 
overlapping subdomains representing the wall, the elastomer and the slab, denoted by 
£2i, £^2 and £^ 3 , respectively. The corresponding domains are chosen as £2] = [0,1] x 
[0,2.8], £2? = [0,1] x [2.8,3] and £^3 = [0,3] x [3,4]. We again used standard linear 
finite elements with 30702 DOFs for these calculations. 

The material parameters E and v will again range from 10— 100 and 0.1 —0.4. 
Since we aim for large numbers of eigenvalues, we perform our simulations for K = 20. 
Fig. 10 shows that we do not only obtain very good convergence for the eigenvalues 
(left) but also for the corresponding eigenfunctions (right). The error curves chosen to 
be represented in Fig. 10 are representative examples for the eigenvalue and eigenfunc¬ 
tion errors in the wall-slab configuration, while the black lines denote the minimum 
and the maximum of the averaged errors over the /i £ E test , respectively. 




Figure 10: Wall-slab configuration with thin elastomer: RB error of eigenvalues (left) 
and eigenfunctions (right) 

The speed-up is similar to the one analyzed in detail in Sect. 5.5. For the wall-slab 
configuration, we show in Table 3 the computation times in the case of K = 20 eigen¬ 
values. As can be seen the computation of the detailed solution takes 14.03 seconds, 
while the computations of the reduced solutions take between 0.14 and 0.24 seconds, 
depending on the basis size N. This results in a speed-up of 100 for N = 50 to 58 if we 
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take N = 300 for an accuracy of 10 7 . 


K 

20 

N 

50 

100 

150 

200 

250 

300 

Detailed solution 

14.03 

Reduced solution 

0.14 

0.16 

0.18 

0.19 

0.22 

0.24 

Reconstruction 

0.027 

0.035 

0.041 

0.049 

0.054 

0.065 

Speed-up 

100 

87 

78 

73 

63 

58 


Table 3: Timings for the detailed solution and the online calculations (reduced solu¬ 
tion including error estimation; reconstruction of the finite element solution from the 
reduced solution) for a slab-wall configuration in seconds 


5.7 Three-dimensional example: First floor building 

Since we aim to apply our results to the modal analysis for vibro-accoustics of lami¬ 
nated timber structures, as they occur in modern timber buildings, we test the perfor¬ 
mance of our method on a three-dimensional geometry representing the first floor of 
a building. Although wooden structures consist of orthotropic materials, we will use 
isotropic material parameters for the ease of computation. 

Usually different materials are used in the construction of a building. In this case, 
we have three different materials for the walls. More precisely, we assume that the outer 
walls are subdomain one, which consists of one material and that the interior walls can 
be divided into two more subdomains, namely ordinary walls and load-bearing walls. 
Fig. 11 depicts our geometry and the corresponding domains. The material parameters 
E and v range from 100 — 1000 and 0.1— 0.4. We perform our simulations for K = 10 
and use standard finite elements with 20994 degrees of freedom. 



Figure 11: Geometry and subdomains: outer walls in red, inner walls load-bearing in 
green, inner ordinary walls in blue 

The first row in Fig. 12 represents the first eigenfunctions for three different param¬ 
eter sets while the second row represents the corresponding fourth eigenfunctions. We 
used the parameter sets fj, 1 = (200,0.1,800,0.3,400,0.2), /. i 2 = (650,0.36,150,0.25,900,0.11) 
and jit 3 = (800,0.3,500,0.1,200,0.4). It can be observed that the eigenfunctions change 
significantly depending on the parameters while still being approximated very well by 
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our method. Fig. 13 shows the error decay for the k-th eigenvalues (left) and eigen¬ 
functions (right), k £ {1,3,5,7,9}, as well as the minimum and maximum averaged 
errors. We again obtain very good convergence. 



Figure 12: Behaviour of the eigenfunctions depending on parameter variations. Top 
row depicts the first eigenfunction and bottom row the fourth eigenfunction. 




Figure 13: First floor configuration: RB error of eigenvalues (left) and eigenfunctions 
(right) 

The speed-up in the three-dimensional setting is even more significant. For the first 
floor of the building, we show in Table 4 the computation times for K = 10 eigenvalues. 
We observe that the computation of the detailed solution takes 31.59 seconds, while the 
computations of the reduced solutions take between 0.084 and 0.142 seconds, depend¬ 
ing on the basis size N. This results in a speed-up of 376 for /V = 50 to 222 if we take 
N = 300 for an accuracy of 10~ 5 . 

6 Conclusion 

In this paper, we have developed a model reduction framework for parameterized ellip¬ 
tic eigenvalue problems and applied it numerically to linear elasticity. We have derived 
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K 

10 

N 

50 

100 

150 

200 

250 

300 

Detailed solution 

31.59 

Reduced solution 

0.084 

0.096 

0.102 

0.111 

0.125 

0.142 

Reconstruction 

0.021 

0.027 

0.031 

0.035 

0.040 

0.046 

Speed-up 

376 

329 

309 

284 

252 

222 


Table 4: Timings for the detailed solution and the online calculations for the first floor 
in seconds 


an asymptotically reliable error estimator for eigenvalues, even for higher multiplic¬ 
ities, that also facilitates an online-offline decomposition. Several eigenvalues have 
been approximated simultaneously by a single reduced space for the variational ap¬ 
proximation of the eigenvalue problem. Altogether we achieve very effective tailored 
greedy strategies for the construction of efficient reduced basis spaces for the simulta¬ 
neous approximation of many eigenvalues. 
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